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Q\ . 

Resonant tunneling through a quantum dot coupled to superconducting reservoirs in the presence 
^ ' of time-dependent external voltage has been studied. A general formula of the current is derived 

, based on the nonequilibrium Green's function technique. Using this formula photon-assisted quasi- 

■ particle transport has been investigated for the quantum dot connected to superconductors. In 

addition, resonant Andreev transport through a strongly correlated quantum dot connected to a 
normal metallic lead and a superconducting lead is studied. 
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I. INTRODUCTION 



Recent advance of nano-technology has stimulated much interest in the study of quantum transport in mesoscopic 
structures. Among the various mesoscopic structures, devices based on quantum dot (QD) have drawn particular 
attention!!! For the quantum dot structures, electron transport is shown to be affected by the confined electrons, in 
the dot. Very recently, novel tunable Kondo effect has been experimentally found in a single electron transistorofl, 
I \ where the level position of the quantum dot and the tunneling rate were controlled by the external gate voltages, via 
'"O ■ Coulomb blockade. 

The influence of time- varying fields on the transport through a quantum dot structure is one of potentially interesting 
areas. Particular applications related to this physics are photo-electron devices, such as single electron pumps, 
turnstiles, and photon detectors (see Ref. [IJfor a review). A time-dependent potential with frequency fl superimposed 
to dc bias potential can induce additional tunneling processes when electrons exchange energy by, absorbing or emitting 
photons of energy O. This kind of tunneling is known as the photon-assisted tunneling (PATp. Experimental study 
on the PAT in the quantum dot devices, based on a configuration of a quantum dot coupled to normal metal leads 



(A), has been reportecuO. 



> 

o 

i In this paper we wish to study tunneling properties in an interacting quantum dot coupled to superconductors in 
the presence of an external ac voltage. For this purpose we first derive a general current formula for this system. 
t-H ' Using this formula we investigate two examples of time-dependent transport, PAT of quasiparticles and that of the 
Kondo- Andreev resonance. When the quantum dot is coupled to a superconductor (S), the resonant tunneling current 
is strongly affected by the singular BCS density of statestt We find that the singularity of BCS density of states plays 
an important role in the photon assisted quasiparticle transport. If the coupling between the leads and the QD is 
sufficiently weak, the subgap transport is suppressed in the case of small dc bias, but the photon assisted tunneling 
combined with the asymmetric BCS density of states could allow finite electron transport. We also investigate 
the many-body resonance by an ac applied voltage in the N-QD-S system as an application of the derived current 
formula. We find new side peaks related to the photon absorption and/or emission of the strongly correlated electrons, 
O ■ m addition to the zero-bias peak due to the Kondo resonance with the Andreev reflection!!! 

^ ' The paper is organized as follows. In Sec. II we derive a time-dependent current formula for the quantum dot 
, connected to superconducting leads. This is an extension of the current formula for the normal metallic leads, derived 

by Jauho, Wingreen, and MehO. On this base, the photon-assisted quasiparticle pumping current is studied in Sec. 

III. The Andreev current for photon-assisted tunneling is studied in Sec. IV. Main results of this paper are summarized 
■ in Sec. V. 



II. CURRENT FORMULA 



We start with a general model Hamiltonian 

ff=V iUcL;Qc4 + E H2 + H D {dl;d a }, (1) 
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where a(s {L, R}, denoting the left hand and right hand sides) is the index of the leads. H a , Hp, and — 
^kcrealVkcrCkcrda + h.c.] represent the Hamiltonians of superconducting (or normal) lead, an interaction region, and 
tunneling between the lead and the interaction region, respectively. The superconducting lead a is characterized by 
the chemical potential [i a and the gap energy A Q . d\ id a ) is the electron creation (annihilation) operator with spin 
a in the quantum dot. When there is a time-dependent voltage difference-fiF" (t) — n a {t) — //d(*) between the lead 
a and the quantum dot, it is convenient to perform gauge transformatior£3. The electron creation and annihilation 
operator for the lead takes the form cl = e*^°/ 2 ct and c ka = e~ l ^ a ^ 2 Ckcr, where d<p a /dt — 2eV a /%. 

The current flowing into the quantum dot can be defined as the rate of change in the number of electrons in a lead 
a. The commutator of the number operator N a = ~^2 kaea ^ ka Cko with the Hamiltonian (|l|) gives rise to the current, 

J a (t) = -!fRe{ £ VZMh*)}. (2) 

In the case of superconducting lead, several kinds of tunneling processes, i.e, quasiparticles, Cooper-pairs, etc., can 
be included in the general current expression. In order to treat the system in a convenient way, we have adopted 
the particle-hole Nambu representations in the lead and the quantum dot such as *& kcr {t) — (c\, a (t) c_fc_ CT (i)) and 
^a{t) = {d\,(t) d-a-it)), respectively. In the particle-hole space, new Green's function can be defined in terms of 
Keldysh Green's function, 



Gj^(*,i , )=i<*i(0® (3a) 

_ /< 4(f) ) ( d- a {t) c ka (t) y 



(3b) 



Ut) = -f Re{ VSG^(i,*)} ii , (4) 



The time-dependent current flowing out of the lead a into the quantum dot can be written as 

2e. 

where {• • -}n stands for the (1,1) component of the current matrix and the hopping matrix is V^ CT = 

Z o \ 

By means of Keldysh techniquefi^l for the noninteraction Hamiltonian of the leads, we have obtained Dyson's 
equation 

- J ^i( g L(i^i)VLG<(ti,i')+g^(i^i)VLG^(ti,t')) , (5b) 

where S — Tc{exp[— i J c drHT(T)]} is the contour-ordered S'-matrix and Tq is the contour-ordering operator. 
G<(t,t') =i( i &i(t')<gi i &Jt)) and G%(t,t') = iQ(t' a (t), *£(*')}) are the lesser and advanced Green's matrices 
in the interaction region. The full Green's functions G a in the interaction region need to be solved. The Green's 
function gfe CT represents the unperturbed Green's function in the lead, which includes the time-dependent phase due 
to the chemical potential difference and the electron level variation. 

In order to obtain the time-dependent Green's functions in the superconducting leads, we take out the time- 
dependent phase associated with a voltage difference in the electron operators as we mentioned above. Assuming A a 
to be independent on the external fields in the superconducting leads, we can separately consider the time variations 
of the energy levels in the leads and in the interaction region (for example, £&(£) = £fc + V a " cos fit ). The total phase 
4> a (t) includes tbe. phase due to the time variation of the energy levels. The Green's functions of a superconducting 
lead have a formliil such as gka(t,t') = R z (<^ a (i))gfco-(i — t')RJ(<^Q,(i' )), where R Z (0 Q ) = cxp[— i((p a /2)er z ] is a rotation 
operator in the two-dimensional complex space. The_electron operators c transform to the quasi-particle operators 7 
by means of the Bogoliubov-Valatin transformationt£l, 



with the BCS coherence factors u k and Vk- The time-dependent Green's functions in the leads can be written as 
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gL(M') = -i9(t-t')RMt)) [e-^^U? +e^(*- t ')V£] R;(<M*'))> (7a) 
g^(t,f) = »R»(^(t)) [/ a (£ fc )e-^*-*')U£ + [1 - /a(^)]e iB ^*-*')Vg] R^i')), (7b) 

where f a {E k ) is the Fermi distribution function in the lead a. One can define the matrix coherence factors U£ and 
V£ given by 

Ug=( < - U f k ) and V£=( ^ V (8) 

The elements of the matrix coherence factors are given by u\ = l/2{l + e k /E k ) and i> 2 = 1/2(1 — e k /E k ). Here, -Bfc = 
\/£ 2 + | A Q | 2 . For normal metallic leads (A Q = 0), £fe becomes the electron energy, and the matrix coherence factors 
are reduced to U£ = 6(e k )P++6(-e k )P- and V£ = 6>(-e fc )P + +6>(e fc )P_, where P± = l/2(l±er 2 ) are the polarization 
matrices and (T z (l) represent the Pauli (unit) matrices. The time-dependent coupling matrix can be written in a 
compact form as T<§(E k ,t,t') = 27riV£(£ fe )V^R z (i)WgR*(i')V£ CT , where W£(E k ) = U£(e fc , E k ) + Vg(-£ fe , -E k ) 
and Ng{E k ) is the density of states of the superconducting lead a. As a consequence, the time-dependent current for 
the systems coupled to two superconducting leads can be written as 

Ja(*) = y f *iy^Im{^e-^( t -* 1 )r§(^,t,t 1 )[-G<(t 1) t) + / a (S fc )G^(t 1 ,t)]} ii , (9) 

with the time-dependent coupling matrix 

S[ ' ' ' ~ ^W~Aj \ e^*)/ 2 ) \ s\gn{E)A a | E | J ^ e- 1 ^')/ 2 ^ 1 Uj 

The coupling constant, which is assumed to be independent of energy and spin, is defined by r a = 2ir \ V k a a | 2 iV"(0) 
with iV Q (0) being the normal density of states at the Fermi level. 

In_the limit A Q = 0, the current expression given in Eq. (|^) reduces to the formula obtained by Jauho, Wingreen and 
MeirEa for normal metallic leads. In principle, the current formula of Eq.([)|) can describe various transport processes, 
associated with the superconducting leads. The time dependent transport physics for a quantum dot connected to the 
superconducting leads with the applied time-dependent voltage can be studied based on Eq. (^|) once the interaction 
in the quantum dot is determined. 



III. PHOTON-ASSISTED QUASIPARTICLE CURRENT 



On the basis of the formula derived in the previous section, we now study the resonant quasiparticle tunneling 
process in the presence of an oscillating external ac voltage. We assume that the quantum dot is weakly coupled 
to the superconducting leads. The Coulomb charging energy (U ~ e 2 /C), where C is the total capacitance of the 
quantum dot, and the level spacing of the discrete electronic states are much greater than the energy gap. Therefore, 
in the weak coupling limit T <C A a -^JJ^Jiie subgap transport due to Andreev reflection is negligible due to large 
charging energy in the quantum dotclllj'EZI. Recently, Oosterkamp, Kouwenhoven, Koolen, Vaart, and HarmansEl 
observed photon induced pumping of the dc current via the 0D state in the normal single electron transistor. They 
also observed photon sideband resonances. Compared to the case of normal metallic leads, here in the case of the 
superconducting leads much enhanced currents are shown in the photon-induced pumping and the direction of the 
current changes, due to the singular BCS density of states. 

We have calculated the time- averaged current, in the presence of sinusoidal external voltage difference between the 
lead a and the quantum dot with the frequency f2 and the oscillation amplitude V" c . By neglecting the off-diagonal 
terms of the general current expression which are related to the Andreev transport, and averaging with respect to 
time, the dc current can be written as 

J S =f E/ ^{rZ(e)g<(e) + K%(e)(g r a (e)-g:(e))}, (11) 

a 

where Q denotes the (1,1) component of the Green's function matrix G. -Tg^e) = J de' r a Qg(e') Q~(e — e') and 
K% (e) = / de' r a g<§ (e') f a (e') Q-(e-e'). The dimensionless BCS factor is g s a (e) = | e \ / ^Je 2 ~ \ A a | 2 . The effects of 
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photon absorption and emission processes are included in the probability function Q^(e) = J^n ^n(^a)S(e ±fi a ±nO) 
where J n {X a ) is the n-th order Bessel function with the argument of the normalized oscillation amplitude A Q = 

ev* e /m. 

To describe the interacting quantum dot we consider the Anderson impurity model Ho — Y] c e a d\d a + Un^ni 
for e a = E- a . In general, the Green's functions of the quantum dot can be found from Dyson's equation Q r ' a — 
£o' a [l + Y, r ' a g r > a ] and G<^. [1 + G r Y, r ]g<[\ + £ Q £ a ] + g r Y><g a . By using the equation of motion technique and the 
mean field approximation!!^, for the photon-assisted tunneling, we obtain the approximate Green's function G„ (e) — 
-{K s {e)/r s {e)){g r a {e) - ^(e)), where r s (e) = E a r s( £ ) and K s(e) = E a K%(e). This approximation guarantees 
automatically the current conservation. The time averaged current through the system becomes 

JQ = JE/ der(s)A$(e), (12) 

where f (e) = (K^{e)r^(e) - Kg(e)rg(e)) /r s (e) and A$(e) = -Im[££(e)] /tt. In the absence of the ac voltage, 
Eq.(|l2[) reduces to the current formula derived in Ref. [I?] for the quasiparticle resonant-tunneling of a quantum dot 
connected to superconductors. After neglecting the higher-order correlation functionsE3, the approximate Green's 
function G^is) can be written as 

1 - (n-a) y f 131 

a{> e- e a -U- E5(e)' K 1 

where the self-energy is £5(e) = — *-Ts(£)/2- Here we do not take into account the Kondo-like correlations because 
it is not important in the case of both leads being superconductors. The occupation number in the nonequilibrium 
state can be obtained from the relation {ri- a ) = (l/2ni) J de g^(s). Then, the current and the occupation number 
of the quantum dot are solved self-consistently. 

When the quantum dot is coupled to the superconductors, the pumping current due to quasiparticles can flow even 
when the ac-fields are applied to both sides of the leads. In the case of the normal lead, o»,tiic other hand, the electron 
pumping current can flow only when the ac-field is applied between the dot and one leadEtEa. In the normal lead case, 
when the ac-fields are applied to both sides of the leads, the current can not flow. Figure [l] displays the zero-bias 
current for unequal superconducting gap energies in the two leads. In this case, since the external frequency is smaller 
than the gap, the single photon processes are suppressed, and two photon processes become the dominant one because 
of the BCS gap. The sharp peak of the current induced by the oscillating field reflects the singularity of the BCS 
density of states. The sharp peaks can give rise to a good spectroscopic resolution. This suggests that a S-QD-S 
system can be utilized as a potentially good photo-electron device. Electron transmission processes are depicted in 
the inset of the Fig.[j]. Note that the direction and the magnitude of the current change as the level position in the 
quantum dot varies. 

Figure || shows photon pumping of the current for two identical superconductors, when a small dc voltage is applied 
to the system. One can clearly see the negative current in some regions of parameters (A, D), which means the 
negative differential conductance arises due to the photon pumping of the current. The dominant tunneling processes 
for Vi = 1.5 A are depicted in the inset for several values of e a . In regions A and B the two-photon process becomes 
dominant, whereas a single photon photon processes dominates, giving rise to a greater current, if the dot level lies in 
between the BCS gaps (C, D). It is noteworthy that the direction and the amplitude of the current change depending 
on the external frequency. As one can see in the Figj^, the direction of the current differs for different frequencies of 
fl = 0.5A, 1.0A and 1.5A. 



IV. PHOTON-ASSISTED ANDREEV CURRENT 



In the present section we will investigate the Andree"v-iy3aonant tunneling through a strongly correlated quantum 
dot coupled to a normal metal and to a superconductorErEl We consider that an ac voltage is applied between the 



normal metallic lead and the quantum dot only. Generalizing the ansatz of Ngfj to the case of a superconducting 
leadEil, (which leads to current conservation,) the time averaged current in the presence of sinusoidal external ac 
voltage on the normal lead can be obtained from Eq.(||) considering the Andreev reflections, 

J£ = jj §Tr{ 0i G;( £ )[S'-(E),K«( e )]G«( £ )}, (14) 

where K^-(e) = [ ^v( £ ) j w jt n K^(e) = J de' T a f a (e') Q^(e — e'). In obtaining the Andreev current wc 
have assumed that the ac voltage difference is applied only on the normal side and that the noninteracting self-energies 
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are in the limit fl, ksT, eV <C A, U. In the absence of the-ac voltage, our equation of ( [ll] ) can be shown to be reduced 
to the current formula derived by Fazio and RaimondiEil. The, Green's functions can be solved by the equation of 
motion technique, taking into account Kondo-like correlations^. In the infinite U limit the Green's function matrix 
has the form 

(£-£ a --£ r (e))G r a (e) = l-(n), (15) 
where £ a = ( £ ° ^ ) and the self-energy can be written as S(e) = S^e) + Ti r s (e), with S^, and Eg denoting 



v 

the self-energy contribution from the normal metal and the superconductor, respectively. The effects of the PAT and 
Kondo correlations are included in H r N u(e) aa23 

yr iVi-V" 7 2/ x n f d£ k r N [l + f N (e k )} 



We can confirm the time reversal symmetry that S^v 22 (e) = — S^ c 11 (— e). In the large A limit (fl, s <C A), ££. can 
be written a£3 

( r s \n{W/A]/2ir -r s /2 \ 

-r s /2 -r s ln[W/A]/2n ) ' (17j 

The off-diagonal terms are directly related to the Andreev transport, whereas the effect of the diagonal term renor- 
malizes the energy level of the quantum dot as £ a ,-a = £<y,-a + r s ln[IK/A]/27r in the large A limit. 

Differential conductance dJ A /dV as a function of the dc bias is shown in Fig.||. The zero-bias peak in the figure 
originates from the Kondo effect associatedpssrith the Andreev reflection. The enhancement of the zero-bias resonance 
by the superconductor is already studiedELrE3. We find that new small peaks appear near the zero-bias peak, due to 
the PAT of the many-body resonance. The|-side peaks are located where the voltage is the multiple of Jl/2, instead 
f2, differently from the case for N-QD-NESEB. This implies that the many-body resonance in the superconducting 
system changes from the normal Kondo resonance because of the Andreev reflection. The Kondo resonance in the 
superconducting system seems to have an effective charge 2e instead of e due to the proximity effect. This proximity 
coupling seems to allow the side peaks at 2eV — nfl with n being an integer. 



V. CONCLUSION 



We have studied resonant transport through an interacting quantum dot coupled to superconductors in the presence 
of external microwave field. We derived a general formula of the current for this system, adopting the nonequilibrium 
Green's function technique. We have seen that the singularity of BCS density of states plays an important as well as 
interesting role in the quasiparticle PAT. Influence of microwave fields on the resonant Andreev transport through a 
strongly correlated quantum dot has also been studied. In addition to a zero-bias peak due to many-body resonance, 
side peaks are found at the bias voltage corresponding to the multiples of the half of the ac oscillation frequency. 
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FIG. 1. Zero-bias quasiparticle current of the S-QD-S in units of eF/Ti for asymmetric gap energy; Al ~ 1.0, Ar = 1.5. 
The other parameters taken here are r = F L = r R = 0.01, k B T = 0.1, U = 20, \ L = \r = 1.0, = 1.0, in units of A L . 
Diagrams of the inset displays the dominant two photon processes for several cases denoted by A, B, C and D. 
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FIG. 2. Photon induced pumping of the quasiparticle current for the S-QD-S in units of er/h for a small dc bias with 
two identical superconductors. = -fi R = eV/2, eV = 0.5, r = r L = r R = 0.01, k B T = 0.1, U = 20, \ L = Xr = 1.0, 
Q — 1.5(solid line), SI = 1.0(dotted line), and Q = 0.5(dash-dotted line) in units of A = Al = Ar. Dominant processes with 
fi = 1.5 are depicted in the diagrams of the inset for several values of the energy level in the QD. 
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FIG. 3. The differential conductance for the N-QD-S structure as a function of dc voltage in the units of 4e 2 //i is plotted at 
the low temperature (k B T = 1 x 10~ 4 ). In the absence of the ac field (sold line), the differential conductance has the zero-bias 
peak due to many-body correlation. For the ac fields of Ajv = 1 with frequency £1 = 0.001 (dotted line) and Q — 0.0015 
(dash-dotted line), the side peaks at 2eV = nQ, are found in addition to the reduced zero- bias peaks. Here r N = T s = 0.02, 
e a = —0.04 on the scale of bandwidth W . 
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